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Abstract 

In  studies  to  characterize  the  kinetic  parameters  for  catalytic  hydrolysis  of  the  organophosphorus 
nerve  agent  1,2,2-trimethylpropyl  methylphosphonofluoridate  (GD)  by  human  paraoxonase 
(HuPONl),  it  was  found  that  standard  analysis  techniques  could  not  be  employed  because  the 
kinetic  parameters  for  the  GD  isomers  were  not  significantly  different  from  each  other. 
Separation  of  the  racemic  mixture  into  component  stereoisomers  for  individual  study  is  feasible 
but  extremely  expensive.  It  was  recognized  that  a  theory  was  required  that  would  allow 
estimates  of  kinetic  parameters  of  each  isomer  to  be  obtained  under  conditions  of  simultaneous 
hydrolysis  of  all  four  stereoisomers.  This  report  presents  a  synopsis  of  the  derivation  of 
theoretical  solutions  for  the  time-course  of  the  four  substrates  (the  stereoisomers)  undergoing 
simultaneous  hydrolysis  by  a  single  active  site  of  an  enzyme.  The  enzyme  “steady  state” 
assumption  made  in  the  derivation  is  tested  by  comparing  the  theoretical  solution  against 
solutions  obtained  by  high  accuracy  numerical  integration  of  the  state  equations.  The  theoretical 
solutions  are  used  to  estimate  the  kinetic  parameters  for  racemic  GD  hydrolysis  by  human 
PONl.  The  graphs  of  the  theoretical  solutions  onto  Lineweaver-Burke,  Eadie-Hofstee,  and 
Hanes-Woolf  plots  provide  visual  demonstrations  of  the  inadequacy  of  those  techniques  for  the 
determination  of  kinetic  parameters  for  multiple  substrate  systems. 


Introduction 

The  reaction  system  to  be  modeled  consisted  of  four  substrates  (the  stereoisomers  of  the  nerve 
agent  soman,  GD)  being  simultaneously  hydrolyzed  by  an  enzyme  (human  paraoxonase)  to 
produce  four  hydrolysis  products  as  depicted  in  Figure  1  shown  below,  which  describes  the 
relationships  among  the  concentrations  of  the  components.  Each  substrate  (A,  B,  C,  D) 
combines  with  free  enzyme  (E)  according  to  a  bimolecular  rate  constant  (kj  kzt,  k7,  kio)  to  form 
an  enzyme-substrate  intermediate  (Ea,  Eb,  Ec,  Ed)  that  progresses  either  forward  according  to  a 
first  order  rate  constant  (k3,  k6,  k9,  ki2)  toward  the  hydrolyzed  substrate  (P,  Q,  R,  S)  and  free 
enzyme  or  backward  according  to  a  first  order  rate  constant  (k2,  k5,  kg,  kn)  to  give  substrate  and 
free  enzyme.  The  enzyme  is  assumed  to  initially  start  entirely  as  free  enzyme  with  concentration 
Eq.  Conservation  of  mass  dictates  that  Eq  =  E  +  Ea  +  Eb  +  Ec  +  Ed  at  any  time. 


Figure  1.  Reaction  schematic  of  the  racemic  GD-HuPON  system.  The  stereoisomers  of  GD  are 
shown  as  A,  B,  C  and  D;  the  free  form  of  HuPON  as  E;  the  GD-PON  intermediates  by  Ea,  Eb, 
Ec  and  ED;  and  the  hydrolysis  products  by  P,  Q,  R  and  S.  Bimolecular  and  first-order  rate 
constants  are  represented  by  the  symbols  ki  through  ki2. 
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Derivation  of  GD-PON  Intermediate  "Steady  State"  Levels 

Mass  balances  on  the  enzyme-substrate  complexes  result  in  the  state  equations  for  the  enzyme- 
substrate  intermediates. 
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After  the  enzyme  “steady  state"  assumption  (that  the  rates  of  change  of  the  intermediates  are 
zero  after  enzyme  steady  state  has  been  achieved),  the  following  relationships  are  obtained: 
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where  the  free  enzyme  level  has  been  written  as  (Eo  -  Ea  -  Eb  -  Ec  -  Ed).  The  equation  set  can 
be  divided  by  Eq  and  rearranged  to  write  it  in  the  form 
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which  may  be  expressed  as  the  matrix  equation 


U  =  M  E  (equation  1) 

4x1  4x4  4x1 


where  U  =  (1,1,1,1)T;E  =  (EA/E0,  EB/E0,  Ec/E0)  ED/E0)T;  and  the  matrix  M  is  given  by 

r  1 

[  1+a  1  1  1  [ 

M  =  j  1  1+b  1  1  j 

[  i  i  i+c  i  j 

[Ill  l+d  | 

L  J 


where  a  =  (k2+k3)/(kiA)  =  KmA/A;  b  =  (k5+k6)/(k4B)  =  KmB/B;  c  =  (k8+k9)/(k7C)  =  KmC/C; 
and  d  =  (kn+ki2)/(kioD)  =  KmD/D. 
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The  inverse  of  M  is 


M  1 
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abed  +  bed  +  acd  +  abd  +  abc 

Pre-multiplying  equation  1  by  M'1  results  in  a  solution  for  the  vector  E  of  steady  state  enzyme 
intermediate  levels  as  functions  of  the  substrate  concentrations  and  the  system  parameters. 

E  =  (bed, acd, abd, abc) T  (abcd+bcd+acd+abd+abc) 

4x1 


Replacing  the  symbols  a,  b,  c,  and  d  using  their  definitions,  we  express  the  fraction  of  enzyme 
existing  in  each  intermediate  form  at  enzyme  steady  state  in  terms  of  the  system  parameters 
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The  fraction  of  free  enzyme  existing  at  enzyme  steady  state  is  given  by 
E/E0  =  1  /(I  +  A/tW  +  B/KmB  +  C/KmC  +  D  /  KmD ) 

Rates  of  Substrate  hydrolysis 

We  write  the  rates  of  change  of  the  substrates  at  enzyme  steady  state 

dA/dt  =  -  ki  A  E  +  k2  Ea 

dB/dt  =  -  k4  B  E  +  k5  Eb 

dC/dt  =  -  k7  C  E  +  k8  Ec 

dD/dt  =  -  ki0D  E  +  knED 


and  substitute  the  steady  state  enzyme  levels.  After  rearrangement  we  find  that 
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(equations  2) 
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The  right  hand  sides  of  the  four  equations  are  the  same,  so  we  may  equate  all  pairs  of  the  left 
hand  sides  to  obtain  6  differential  relationships  between  the  substrates 
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These  relationships  can  be  integrated  with  respect  to  time  from  0  ->  T  to  obtain  relationships 
between  the  substrates  at  any  time 
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And  where  VmaxA  —  kyEo,  Vmaxg  —  kgEo,  Vmaxc  —  kyEo,  VmaXD  —  k  i  oEq 


Using  these  relationships,  we  replace  the  variables  in  equations  2  in  favor  of  each  equation’s 
substrate  and  rearrange  the  resulting  equations  to  obtain 
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Where 

P  —  (KmB/KmA)  (  Vmax^/ VmaxB  )  }  cj  —  (Kmc/Km^)  ( Vmax^/ Vmaxc )  }  T  —  (K  mo/KmA.)  (  Vmax^/ Vmaxg  ) 

s  =  ( KmC / KmB )  (VmaxB/VmaxC)  ;  t  =  (KmD/KmB)  ( VmaxB/VmaxD)  ;  u  =  (K  mo/Kmc)  ( Tna  xC  /  Ana  xD  ) 


(equations  3) 


The  above  equations  each  contain  a  single  variable  and  are  of  the  form 

constant  =  F  {  x(t)  }  d/dt{  x(t)  }, 

which  can  be  integrated  with  respect  to  time  to  obtain  functions  of  the  form,  }F{x}dx  = 
constant*T.  Each  integral  can  be  evaluated  on  the  interval  0  to  T  to  obtain  equations  relating 
time  to  the  fractional  substrate  at  that  time  level. 
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For  the  first  equation  in  equation  set  3, 


Vmaxfl/KmA  —  (  (A/Aq)  1  +  A()KmA  +  B0  (A/A0)  <1/p)-1  tW1 
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Integrate  each  term  with  respect  to  time  from 


'Vmaxfl/KmA  dt  =  /  (A/A0)  M/dttA/Aol  dt  + 

S  AoKmA 

d/dt  {A/ A0 } 

dt 

+ 

f  B0K ms-1 

(A/A0) 

<1/p)_1  d/ dt  { A/A0 } 

dt 

+ 

f  Coiwr1 

(A/A0) 

(1/q)_1  d/dt{A/A0} 

dt 

+ 

f  DoV1 

(A/A0) 

<1/r)~1  d/dt{A/A0} 

dt 

And  evaluate  at  t=T  and  t=0 


t=T  t=T  t=T 


-VmaxA/KmA  t]  |  =  [  Ln  { A<t  >  /  A0 }  ]  1  + 

t=0  t=0 

[  AqK^-1 

A<t>/A0 ] 

t=0 

t=T 

+ 

[  BoKmB 

(p)  (A<t>/A0)  <1/p>  ] 

t=0 

t=T 

+ 

[  C0Km(f 1 

(q)  (A<t>/A0)  (1/q)  . 

]  1 

t=0 

t=T 

+ 

[  BoKmD 

(r)  (A<t>/A0)  (1/r)  . 

]  1 

t=0 

to  obtain 


"Vmaxfl/Kjnft  T  -  0  =  Ln{A<T>/A0}  -  Ln{A<0>/A0 

+  AgK^-1  A<T>/A0  -  AqK^-1  A<0>/Ag 
+  p  B0KmB_1  (A<T>/A0)  <1/p)  -  p  B0KmB_1  (A<0>/A0)  <1/p) 

+  q  CoKmC_1  (A<T>/A0)  <1/q)  -  q  CoK^1  (A<0>/A0)  (1/q) 

+  r  D0K„Dq  (A<T>/A0 )  <1/r)  -  r  DqPW1  (A<0>/A0)  (1/r) 

— WaxA/KmA  T  =  Ln{A/Ao}  +  AqK^  A/ Aq  AoK^ 

+  pB0KmB  1  (A/A0)  <1/p)  -  B0KmB  1 
+  qCoKmC_1  (A/A0)  <1/q)  -  CqK^1 
+  rD0KmD_1  (A/A0)  <1/r>  -  D0KmD_1 
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We  replace  p,  q,  and  r  using  their  definitions  (but  not  in  the  exponents)  to  write 


WaxA/KmA  T  -  Ln{A/A0}  + 

AofW1  { 

A/A0  -  1  } 

+ 

(WIU) 

(^maxA/^maxB )  ^0-^-mB 

{ 

(A/A0)  <1/p) 

-  1 

} 

+ 

(Kmc/Kmft) 

( ^maxA  /  VmaxC  )  C  0  Kmc 

{ 

(A/A0)  <1/q) 

-  1 

} 

+ 

( KmD  /  Km* ) 

( ^maxA  /  VmaxD  )  D  q  Kmpj 

{ 

(A/A0)  <1/r) 

-  1 

} 

T  =  -(WV^ILnlA/Aol 

_  7\  JZ  “1 


AoKmA 

{ 

A/A0 

1 

—  (KjnA/VmaxA) 

(Kme/Km*) 

(^maxA/^maxB )  -^O^mB 

{ 

(A/A0)  <1/p)- 

1 

( ^-mA  /  ^maxA ) 

( KmC  /  KmA ) 

("^maxA/^maxC )  ^O^mC 

{ 

(A/A0)  (1/q)  - 

1 

( ^mA  /  ^maxA ) 

( KmD  /  KmA ) 

(^maxA/^maxo)  ^O^mD 

{ 

(A/A0)  (1/r)- 

1 

From  which  we  obtain 


Ta 

( Ao/VmaXA) 

(1  - 

(A/A0) 

(KmA/KmA)  (VmaxA/VmaxA)  ^ 

+ 

( Bo  /  VmaxB ) 

(1  - 

(A/A0) 

( Km*/ KmB )  (VmaxB /VmaxA)  ^ 

+ 

(C0/Vmaxc) 

(1  - 

(A/A0) 

(Kma/Kmc)  (VmaxC/VmaxA)  j 

+ 

(Do/VmaxD) 

(1  - 

(A/A0) 

(Km*/ KmD)  (VmaxD/VmaxA)  ^ 

— 

( KmA  /  WnaxA ) 

Ln  (A/A0) 

Similarly,  we  obtain  from  the  other  three  equations 

Tb 

+ 

( Ao/VmaxA) 

(  1  - 

(B/B0 

^  (KmB/KmA)  (VmaxA/VmaxB) 

+ 

(Bo/VmaxB) 

(  1  - 

(B/B0 

^  (KmB /KmB)  (VmaxB/VmaxB) 

+ 

(Co/VmaxC) 

(  1  - 

(B/B0 

J  (KmB/Kmc)  (VmaxC/VmaxB) 

+ 

(Do/VmaxD) 

(  1  - 

(B/B0 

)  (KmB /KmD)  (VmaxD/VmaxB) 

- 

( KmB  /  VmaxB ) 

Ln  (B/B0) 

Tc 

( Ao/VmaXA) 

(  1  - 

(C/Co 

^  (Kmc/KmA)  (VmaxA/VmaxC) 

+ 

(Bo/VmaxB) 

(  1  - 

(C/C0 

J  (Kmc/Kms)  (VmaxB/VmaxC) 

+ 

(C0/Vmaxc) 

(  1  - 

(C/C0 

J  (Kmc/Kmc)  (VmaxC/VmaxC) 

+ 

(Do/VmaxD) 

(  1  - 

(C/Co 

J  (Kmc/KmD)  (VmaxD/VmaxC) 

- 

(Kmc/Vmaxc) 

Ln  (C/C0) 

Td 

( Ao/VmaXA) 

(  1  - 

(D/D0 

)  ( KmD / KmA )  (VmaxA/VmaxD ) 

+ 

(Bo/VmaxB) 

(  1  - 

(D/D0 

J  ( KmD / KmB )  (VmaxB/VmaxD) 

+ 

(C0/Vmaxc) 

(  1  - 

(D/D0 

)  (KmD/Kmc)  (VmaxC/VmaxD ) 

+ 

(Do/VmaxD) 

(  1  - 

(D/D0 

^  ( KmD / KmD )  (VmaxD/VmaxD) 

- 

( KmB  /  VmaxB ) 

Ln  (D/D0) 

(equations  4) 

Although  complex,  these  equations  give  the  time  it  would  take  for  each  substrate  (normalized  by 
its  initial  level)  to  reach  a  particular  level.  As  such,  they  can  be  used  to  form  graphs  of  the 
substrate  levels  as  functions  of  time  by  treating  the  fractional  substrate  level  (X/X o)  as  the 
independent  variable.  Since  the  enzyme  intermediate  levels  start  at  zero,  the  fractional  substrate 
levels  can  only  fall  during  hydrolysis.  Therefore,  they  each  start  at  a  value  of  one  and  fall 
towards  a  value  of  zero  (most  likely  at  different  rates)  as  the  substrates  are  hydrolyzed. 
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Application  of  the  Theory 

In  studies  where  racemic  GD  was  hydrolyzed  by  human  PON1,  the  activity  of  each  GD  isomer 
was  measured  at  various  times  during  the  hydrolysis  (as  detailed  in  “Direct  detection  of 
stereo  specific  soman  (GD)  hydrolysis  by  human  serum  paraoxonase  [HuPONl]”  by  Yeung  et  al., 
FEBS  Journal,  in  press).  The  functions  giving  the  theoretical  solutions  for  the  four  isomers 
(equation  set  4)  were  plotted  on  an  X-Y  chart  with  the  fractional  substrate  levels  on  the  ordinate 
and  the  time  taken  for  the  substrate  to  fall  to  a  particular  substrate  level  on  the  abscissa. 

Measured  stereoisomer  concentrations  and  times  were  also  plotted  on  the  chart.  The  kinetic 
parameters  of  the  equations  were  adjusted  to  fit  the  theoretical  curves  to  the  actual  data.  The 
SigmaPlot®  and  Prism®  nonlinear  regression  curve  fitting  programs  were  used  in  an  attempt  to 
automate  the  curve  fitting  process,  but  the  methods  could  not  be  made  to  converge,  presumably 
because  of  the  large  number  of  parameters  (12)  and  the  complexity  of  the  functions.  Consequently, 
the  curves  were  fit  manually  by  adjusting  parameters  in  an  Excel®  workbook  that  had  been 
programmed  to  compute  and  chart  the  theoretical  solutions  along  with  measured  data  points. 

A  standard  assumption  when  dealing  with  enzyme-substrate  systems  is  that  the  maximum  rate  of 
product  generation  (and  substrate  consumption)  occurs  at  the  beginning  of  the  enzyme  action. 
While  this  is  true  in  a  single  substrate  system,  in  multiple  substrate  systems  the  hydrolysis  of  a 
particular  isomer  can  start  slowly  because  the  enzyme  is  busy  hydrolyzing  another  substrate. 

The  hydrolysis  of  the  substrate  can  speed  up  when  the  enzyme  becomes  available.  With  this 
concept  in  mind,  the  Vmax  value  for  an  isomer  was  set  to  the  maximum  hydrolysis  velocity 
measured  at  any  time  during  the  hydrolysis  experiment.  Given  these  Vmax  values  and  the  known 
starting  enzyme  concentration,  Eo,  estimates  of  the  kinetic  parameters  (k3,  k6,  k9,  k  1 2)  were  found. 
Under  the  assumption  that  the  single  enzyme  intermediate  could  disassociate  in  the  forward  or 
reverse  directions  with  equal  facility,  the  parameters  (k2,  k.5,  ks,  kn)  were  set  equal  to  the  values 
for  (k3,  k6,  kg,  kn).  This  left  just  the  bimolecular  rate  constants  (ki  k4,  k7,  kio)  to  be  adjusted 
using  curve  fitting.  The  results  for  a  particular  initial  GD  level  are  shown  below  in  Figure  2. 

Hydrolysis  of  0.37  mM  racemic  GD  by  Human  PON1 


0  50  100  150  200  250 

Time  (minutes) 

Figure  2.  Curve  Fitting  of  the  Theoretical  Solutions  to  Measured  Data 
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Comparison  with  Numerical  Simulations  -  Does  the  theory  work? 

To  test  the  theoretical  results,  AcslXtreme  (a  numerical  simulation  software  package)  was  used 
to  provide  high  accuracy  numerical  simulations  of  the  original  system.  The  results  of  the 
simulations  were  plotted  on  the  X-Y  Chart  shown  in  Figure  3. 


Comparison  of  Theoretical 
and  Numerical  Solutions 


—  C(-)P(-)  theoretical 
C(-)P(-)  numerical 

—  C(-)P(+)  theoretical 
C(-)P(+)  numerical 
C(+)P(-)  theoretical 

—  C(+)P(-)  numerical 

—  C(+)P(+)  theoretical 
C(+)P(+)  numerical 


100  150 

Time  (minutes) 


300 


Figure  3.  The  theoretical  solutions  for  the  four  substrates  (drawn  in  bold)  are 
shown  with  the  numerical  solutions.  The  solutions  are  identical. 


To  explore  the  amount  of  error  introduced  by  the  enzyme  steady  state  assumption,  the  numerical 
levels  of  the  enzyme  intermediates  were  plotted  along  with  the  theoretical  values  for  them 
predicted  from  the  enzyme  steady  state  assumption.  The  chart  is  shown  in  Figure  4. 
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Comparison  of  Assumed  Enzyme  "Steady  State" 
Levels  and  Actual  (Numerically  Integrated)  Levels 


0.00001  0.0001  0.001  0.01  0.1  1  10  100  1000 

Time  (minutes) 


Figure  4.  The  effect  of  the  enzyme  steady  state  assumption.  Note  that  the  differences 
between  the  assumed  and  actual  levels  vanish  quickly. 


What  about  the  single-substrate  techniques? 

It  is  illustrative  to  examine  what  the  theoretical  solutions  would  look  like  on  Lineweaver-Burke, 
Eadie-Hofstee,  and  Hanes-Woolf  plots.  To  accomplish  this,  we  used  the  negative  reciprocal  of 
the  rate  of  change  of  a  substrate  as  the  “1/Velocity”  variable.  The  negative  is  used  because  the 
rate  of  the  hydrolysis  is  the  negative  of  the  rate  that  the  substrate  changes.  The  reciprocals  of  the 
equations  describing  the  rates  of  change  of  the  substrates  were  rearranged  and  the  relationships 
between  the  substrates  were  used  to  write  each  equation  in  terms  of  only  a  single  substrate. 


1/VA 


1/Vb 


1/Vc 


1/V; 


maxA 


l/V; 


maxB 


1  /Vmaxc 


+  (WVmaxA)  (1/A) 


+ 

( Vmax^Ao ) 

(B0/KmB) 

(A/A0) 

(KmA^maxB/  (KmBVmaxA) 

-  i) 

+ 

(  Kjha/  VmaxA-^-0  ) 

(C0/Kmc) 

(A/Aq) 

(KmAVmaxc/  (KmcVmaxA) 

-  i) 

+ 

(  Kjtia/ VmaxA-^-0 ) 

(D0/KmD) 

(A/Aq) 

(KmA^maxD/  (KmDVmaxA) 

-  i) 

+ 

(KmB/VmaXB)  (1/B) 

+ 

( KmB  /  V maxB^  o  ) 

(Ao/fW) 

(B/B0) 

(KmB^maxA/  ( KmA^maxE ) 

-  i) 

+ 

(  KmB  /  V maxBB  o  ) 

(C0/Kmc) 

(B/B0) 

(KmBVmaxc/  (KmcVmaxB) 

-  i) 

+ 

( KmB  /  V maxBB  o  ) 

(D0/KmD) 

(B/B0) 

(KmBVmaXD/  (KmDVmaXB) 

-  i) 

+ 

(KmC/Vmaxc)  (1/C) 

+ 

(KmC/VmaxcCo) 

(Ao/IW) 

(C/C0) 

(KmC^maxA/  ( KmA^maxC ) 

-  i) 

+ 

(KmC/VmaxcCo) 

(B0/KmB) 

(C/C0) 

(KmC^maxB/  (KmBVmaxc) 

-  i) 
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1/VD 


1/V; 


maxD 


+  (KmC/VmaxCCo)  (D0/KmD)  (C/Co)  <KmcVmaxD/  (KmDvmaxC)  -  U 
+  (KmD/VmaxD)  (1/D) 

+  (KmD/VmaxDDo)  (Ao/Kn*)  (D/Do)  <KraDVmaxA/ (W^xd)  -  D 
+  (KmD/VmaxDD0)  (B0/KmB)  (D/D0)  (k^v»-b/ "  2) 
+  (KmD/VmaxDD0)  (C0/Kmc)  (D/D0)  <KmCvmaxD)  -  D 


(equations  5) 

The  first  two  terms  on  the  right  of  each  equation  give  the  standard  (linear)  L-B  equation  as  it 
would  apply  to  a  single  substrate.  The  last  three  terms  are  power  functions  of  the  fractional 
substrate  levels.  The  terms  add  non-linear  contributions  to  the  (1/V)  dimension.  Figure  5 
presents  some  measured  data  and  the  corresponding  theoretical  solutions  plotted  on  the  double 
reciprocal  plot. 


Lineweaver-Burke  Plot  of  Theoretical  Solutions 
and  Measured  Data  for  Hydrolysis  of 
1.67  mM  Racemic  GD  by  Human  PON1 
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Figure  5.  The  Lineweaver-Burke  plot  for  theoretical  solutions  of  multi-substrate 
hydrolysis  plotted  along  with  measured  data  points.  Note  the  nonlinearity  for  the 
higher  substrate  concentrations. 


On  this  chart,  each  isomer  starts  at  its  highest  (initial)  concentration  that  occurs  at  the  leftmost 
end  of  its  curve.  As  the  isomer  level  falls,  the  curve  is  traversed  from  left  to  right.  It  can  be  seen 
from  the  figure  that  the  nonlinearity  is  most  pronounced  when  the  isomer  concentrations  are 
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highest.  Note  that  the  reaction  velocities  of  three  of  the  four  isomers  speed  up  (falling  to  lower 
1/V  levels)  initially  as  the  bulk  of  the  enzyme,  which  was  initially  complexed  with  the  fourth 
isomer,  becomes  available.  Since  the  L-B  transformation  does  not  linearize  the  system,  it  is  not 
appropriate  to  use  the  method  to  estimate  the  kinetics  of  a  multiple  substrate  system. 
Multiplying  each  equation  in  the  equations  5  set  by  the  corresponding  substrate  level  gives  the 
theoretical  solutions  in  the  Hanes-Woolf  format. 


a/va  = 

A/  VmaxA 

+ 

( Vmax£.) 

+ 

( KjnA./ Vmax^) 

(B0/KmB) 

(A/A0) 

(KmAVmaxB)  /  (KmBVmaxA) 

+ 

(  Kjha/  Vmax^) 

(C0/Kmc) 

(A/Aq) 

(KmAVmaxC)  /  (KmC VmaxA) 

+ 

( K-mA./ VmaxA) 

(D0/KmD) 

(A/Aq) 

(KmA^maxo)  /  (KmD VmaxA) 

B/Vb  = 

B  /  Vmaxg 

+ 

(  KmB  /  V  maxB ) 

+ 

(  KmB  /  V  maxB  ) 

(Ao/Kha) 

(B/B0) 

( KmB VmaxA)  /  (KmAVmaxB) 

+ 

(  KmB  /  V maxB  ) 

(C0/Kmc) 

(B/B0) 

(KmBVmaxC)  /  (KmcVmaxB) 

+ 

(  KmB  /  V maxB  ) 

(D0/KmD) 

(B/B0) 

(KmBVmaxD)  /  (KmDVmaxB) 

C/Vc  = 

c/vmaxC 

+ 

(Kmc/Vmaxc) 

+ 

(Kmc/Vmaxc) 

(Ao/K„a) 

(C/Co) 

(KmC VmaxA)  /  (KmAVmaxC) 

+ 

(Kmc/Vmaxc) 

(B0/KmB) 

(C/Co) 

(KmcVmaxB)  /  (KmBVmaxC) 

+ 

(Kmc/Vmaxc) 

(D0/KmD) 

(C/Co) 

(KmcVmaxD)  /  (KmDVmaxC) 

D/Vd  = 

D  /  VmaXD 

+ 

(  KmD  /  V maxD  ) 

+ 

(  KmD  /  V maxD  ) 

(  Aq/Kha) 

(D/D0) 

(KmD VmaxA)  /  (KmAVmaxD) 

+ 

( KmD  /  V  maxD  ) 

(B0/KmB) 

(D/D0) 

(KmDVmaxB)  /  (KmBVmaxD) 

+ 

(  KmD  /  V  maxD  ) 

(C0/Kmc) 

(D/D0) 

(KmDVmaxC)  /  (KmcVmaxD) 

(equations  6) 


The  first  two  terms  in  each  equation  give  the  single  substrate  form  of  the  equations  while  the  last 
three  terms  add  nonlinearities  to  the  S/Vs  dimension. 
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Hanes-Woolf  Plot  of  Theoretical  Solutions 
and  Measured  Data  for  Hydrolysis  of 
1.67  mM  Racemic  GD  by  Human  PON1 


(mM) 

Figure  6.  Hanes-Woolf  plot  of  the  multiple-substrate  theoretical  solutions  and  some 
measured  data  points. 


Each  equation  in  equations  5  may  be  multiplied  by  VsVmaxs  and  rearranged  to  form  relationships 
that  can  be  plotted  on  an  Eadie-Hofstee  plot. 


VA 

—  ^maxA  ~~ 

KmA 

{ 

l 

+ 

(B0/KmB) 

(p^/p^  'j  KmAVmaxB/KmBVmaxA 

+ 

(C0/Kmc) 

(p^/ p\  J  KmAVmaxC/KmCVmaxA 

+ 

(D0/KmD) 

(p^/ p^  J  KmAVmaxD/KmDVmaxA 

}  (Va/A) 

VB 

—  VmaxB  — 

K.B 

{ 

l 

+ 

(Ao/fW) 

(~Q  /  ^  KmBVmaxA/KmAVmaxB 

+ 

(C0/Kmc) 

^gy»g  ^  KmBVmaxC/KmCVmaxB 

+ 

( D  o  /  Kmn ) 

^gy»g  j  KmBVmaxD/KmDVmaxB 

}  (Vb/B) 

Vc 

—  ^maxC 

Kmc 

{ 

l 

+ 

(Ao/IU) 

((2/(2  'j  KmCVmaxA/KmAVmaxC 

+ 

(B0/KmB) 

((2/(2  )  KmCVmaxB/KmBVmaxC 

+ 

(D0/KmD) 

((2/(2  )  KmCVmaxD/KmDVmaxC 

}  (Vc/C) 

vD 

—  ^maxD  — 

K-mD 

{ 

l 

+ 

(Aq/K^) 

(Y)  /  Y)  J  KmDVmaxA/KmAVmaxD 

+ 

(Bo/K^) 

^g/g  ^  KmDVmaxB/KmBVmaxD 

+ 

(C0/Kmc) 

^g^g  j  KmDVmaxC/KmCVmaxD 

}  (V„/D) 

(equations  7) 
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Eadie-Hofstee  Plot  of  Theoretical  Solutions 
and  Measured  Data  for  Hydrolysis  of 
1.67  mM  Racemic  GD  by  Human  PON1 


Reaction  Velocity  /  Substrate  Concentration 
(1/min) 


Figure  7.  The  data  and  theoretical  solutions  used  in  the  previous  two  figures  plotted  on 
the  Eadie-Hofstee  plot.  Single  substrate  systems  plot  as  a  straight  line  with 
slope  =  -Km  and  y-intercept  =Vmax. 


Exploring  the  form  of  equations  7  provides  an  insight  into  the  effect  that  substrate  multiplicity 
has  on  the  apparent  affinity  (Km).  The  quantities  in  braces  "  { } "  in  the  equations  collapse  to  the 
value  1  when  the  initial  levels  of  competing  substrates  are  set  to  zero.  The  equations  become  the 
Eadie-Hofstee  equations  for  independent  single  substrate  systems.  The  effect  of  substrate 
multiplicity  is  to  increase  the  apparent  affinity  that  would  be  measured  for  a  single  substrate. 
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Conclusions 

We  conclude  that  the  equations  presented  as  equation  set  4  are  the  correct  theoretical  solutions 
for  the  enzyme-substrate  system  presented  as  Figure  1  based  upon  the  derivations  and  a 
comparison  of  numerical  and  theoretical  solutions.  The  enzyme  steady  state  assumption  made  in 
the  derivation  of  the  solutions  has  little,  if  any,  effect  on  the  solutions  after  an  initial  period  less 
than  one  second.  We  have  found  that  use  of  the  theoretical  solutions  to  estimate  the  kinetic 
parameters  of  the  system  is  possible  and  yields  sensible  results.  We  can  not  conclude  that  the 
system  presented  as  Figure  1  is  an  accurate  representation  of  the  GD-PON  hydrolytic  system 
although  it  does  represent  the  simplest  GD-PON  system.  The  equations  presented  here  are  for  a 
four  substrate  system,  but  systems  of  two  and  three  substrates  have  also  been  solved  and  have 
solutions  that  are  subsets  of  the  presented  work.  A  formal  approach  can  be  used  to  extend  the 
equations  to  any  number  of  multiple  substrates.  In  general,  this  work  should  be  useful  for 
estimating  the  kinetic  parameters  of  any  multiple-substrate,  single  active  site  system.  By  setting 
kinetic  parameters  judiciously,  this  work  should  also  prove  useful  in  the  study  of  enzyme 
inhibitors  (organophosphates)  and  protectants  (carbamates)  that  can  both  be  viewed  as  special 
types  of  competitive  substrates. 
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